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ABSTRACT 


V 

This  paper  considers  laminated  beams  with  layers  of  different  isotropic 
or  orthotropic  materials  fastened  together  by  thin  adhesives.  The 
stresses  that  result  from  subjecting  each  component  layer  of  the  beam 
to  different  temperature  or  moisture  stimuli  which  may  also  vary  along 
the  length  of  the  beam,  are  calculated.  Two-dimensional  elasticity 
theory  is  used  so  that  a  wide  range  of  problems,  such  as  that  of  beams 
composed  of  layers  of  orthotropic  materials  like  wood,  can  be  studied, 
and  accurate  distributions  of  normal  and  shear  stresses  obtained.  The 
stress  intensity  along  the  bearing  surfaces  of  the  layers  of  the  beam 
is  of  particular  importance  because  it  is  responsible  for  delamination 
failures  of  laminated  structural  elements.  The  distributions  of  inter¬ 
laminar  normal  and  shear  stresses  measured  along  the  longitudinal  axis 
of  the  beam  indicate  that  high  stress  intensity  occurs  in  the  end  zones 
of  the  beam.  Thus,  delamination  failure,  when  it  occurs,  will  start  at 
the  ends  of  the  beam. 
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LAMINATED  BEAMS  OF  ISOTROPIC  OR  ORTHOTROPIC  MATERIALS 


SUBJECTED  TO  TEMPERATURE  CHANGE 


By 

SHUN  CHENG-/ 
and 

T.  GERHARDT£/ 


INTRODUCTION 


Beams  are  among  the  most  widely  used  structural  elements.  The  appli¬ 
cation  of  laminated  beams  is  expanding  as  they  are  studied  and  developed 
by  the  forest  products  industry  as  well  as  other  industries  [1,2,3]. 
However,  the  delamination  failure  or  deformation  (bow,  twist,  warp,  etc.) 
induced  in  laminated  beams  by  thermal  or  moisture  stimuli  has  always 
been  of  major  concern. 

The  use  of  elementary  beam  theory  in  solution  of  this  problem  does  not 
allow  evaluation  of  the  shearing  and  normal  stresses  along  the  bearing 
surface.  Thus,  these  stresses  cannot  be  determined  from  Timoshenko's 
pioneering  analysis  [4]  of  a  bimetal  strip  submitted  to  uniform  heating 
along  its  length.  Notable  contributions  on  this  subject  due  to  Boley 
and  others  may  be  found  in  [5,6].  However,  Grimado's  analysis  [7]  is 
a  further  consideration,  extension,  and  significant  improvement  of  the 
same  problem  treated  by  Timoshenko.  When  Grimado  [7]  uses  elementary 
beam  theory  the  effect  of  the  bonding  material  between  the  two  layers 
of  the  strip  is  taken  into  account  by  treating  the  bonding  material  as 
a  third  layer.  Grimado  deduces  a  sixth-order  governing  differential 
equation  [7]  compared  to  the  fourth-order  biharmonic  equation  of  plane 
stress  problems.  It  is  shown  in  [7]  that  the  sixth-order  equation 
reduces  to  a  characteristic  cubic  equation  which,  unlike  the  biharmonic 
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equation,  may  yield  complex  roots  when  solved.  Thus,  the  analysis 
presented  in  [7]  is  not  necessarily  as  simple  and  direct  as  that  based 
on  the  two-dimensional  elasticity  theory  used  here. 

In  the  current  paper  a  beam  of  uniform  cross-section  made  of  layers 
of  different  isotropic  or  orthotropic  materials  (such  as  wood,  fastened 
together  by  thin  adhesives)  is  considered  in  accordance  with  two- 
dimensional  elastic!  y  theory.  Instead  of  being  subjected  to  a  uniform 
heating,  treated  previously  (4,7],  each  layer  may  have  different  tem¬ 
perature  distributions  along  its  length.  The  use  of  two-dimensional 
elasticity  theory  should  accurately  yield  the  distribution  of  shear 
and  normal  stresses  in  the  beam.  The  interlaminar  stresses  between 
layers  are  known  to  be  mainly  responsible  for  delamination  failures 
of  laminated  beams.  If  the  stress  distribution  between  layers  can  be 
determined,  such  failures  may  be  minimized  or  eliminated  by  an 
appropriate  choice  of  materials  and  beam-section  properties.  The  use 
of  two-dimensional  elasticity  theory  can  be  further  justified  in  the 
analysis  of  beams  of  orthotropic  materials  since  in  such  a  case 
elementary  beam  theory  cannot  take  into  account  the  effects  of  material 
properties  on  stresses  and  deformation. 

It  is  difficult  to  accurately  estimate  the  effect  of  the  bonding  material 
on  the  stress  distribution  and  deformation  of  the  laminated  structural 
element  without  also  treating  the  thin  layers  of  the  bonding  material 
in  the  same  way  that  the  component  layers  of  primary  concern  are  treated. 
When  the  bonding  material  between  two  component  layers  is  very  thin, 
the  effect  of  the  bonding  material  on  the  stress  distribution  and 
deformation  of  laminated  beams  can  be  negligible.  Thus,  this  paper 
considers  two  cases:  (1)  beams  having  two  component  layers  rigidly 
bonded  together  and  (2)  beams  having  three  layers  where  the  bonding 
material  is  treated  as  a  third  layer.  These  are  analyzed  so  that  the 
effect  of  the  bonding  material  can  be  revealed.  Special  attention  is 
given  to  the  stress  intensity  along  the  bearing  surfaces  of  the  layers 
of  the  laminated  beams  since  it  is  responsible  for  delamination  failures. 
Only  laminated  beams  subjected  to  thermal  stimulus  are  treated  in  detail. 
However,  beams  subjected  to  moisture  stimulus  can  be  analyzed  in  the 
same  way  simply  by  replacing  the  linear  coefficient  of  thermal  expan¬ 
sion  multiplied  by  the  change  of  temperature  with  a  linear  coefficient 
of  moisture  shrinkage  multiplied  by  a  relative  moisture  content 
[8,  p.  3-11]. 


TWO-DIMENSIONAL  THERMOELASTIC  FORMULATION 


Consider  a  beam  of  unit  width  which  is  made  of  two  layers  of  different 
materials  fastened  together  by  a  thin  adhesive,  where  each  layer  is 
subjected  to  different  arbitrary  temperature  distribution  T^Cx), 
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III  up  I 


(i  =  1,2)  along  its  length.  The  laminated  structure  is  initially  free 
from  stress.  Let  the  x-axis  be  the  longitudinal  axis  which  lies  along 
the  bonding  line  of  the  two  layers  of  the  beam  of  length,  SL,  total 
thickness,  h,  and  place  the  y-axis  at  the  left  end  of  the  beam.  Let 

hj  be  the  thickness  of  the  top  layer,  h^  the  thickness  of  the  bottom 

layer,  and  d}  =  h^,  d^  =  -h^.  Let  E  be  the  Young's  modulus,  G  the 

shear  modulus,  and  v  (Greek  letter  nu)  the  Poisson's  ratio.  Let  a  and 

e  be  the  components  of  stress  and  strain,  u  and  w  the  displacements 

along  the  x  and  y  directions,  and  a  and  a  the  coefficients  of 

x  y 

thermal  expansion  in  directions  x  and  y,  respectively.  The  next 
section,  concerning  laminated  beams  of  isotropic  materials,  is  limited 
to  the  consideration  of  a  specific  case  where  the  beam  is  uniformly 
heated  by  raising  its  temperature  t  degrees.  It  is  a  fundamental  case 
that  can  reveal  the  essential  features  of  thermal-stress  problems  of 
laminated  beams;  the  solution  can  be  used  to  make  comparisons  with  known 
results  in  the  literature.  A  subsequent  section,  concerning  laminated 
beams  of  orthotropic  materials,  treats  the  general  case  of  temperature 
distribution . 


LAMINATED  BEAMS  OF  ISOTROPIC  MATERIALS 


In  this  section  we  consider  a  laminated  beam  consisting  of  two  different 
isotropic  materials.  For  plane  stress  distribution  the  stress-strain 
and  strain-displacement  relations,  the  governing  differential  equation, 
and  the  expressions  for  stresses  in  terms  of  a  stress  function  are  [9,10]. 
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When  the  laminated  beam  is  uniformly  heated,  changing  the  temperature 
by  t  degrees,  may  be  expressed  as 


T. 

1 


os 

z 

n 


1 

n 


sin 


nnx 

£ 


(n  =  1,3,5,...) 


(4) 


Accordingly,  the  biharmonic  function  <(>^  for  the  ith  layer  which 
satisfies  equation  (2)  may  be  taken  as 


os 

d>.  =  Z  (A  .cosh  vy  +  B  ,vy  sinh  Yy  +  C  .sinh  Yy 
l  ,  ni  m"  13  ni  13 

n=l 

+  Dniyy  cosh  yy)sin  yx  (n  =  1,3,5...)  (5) 


where  Y  =  — . 


(6) 


The  boundary  conditions  on  the  top  and  bottom  surfaces  of  the  beam  are 

a  =  0,  t  =0,  on  y  =  d.(d,  =  h, ,  d  =  -h„)  (7) 

y  xy’  3  i  1  12  2 

In  addition  to  the  preceding  boundary  conditions,  the  conditions  of 
continuity  along  the  line  of  division  of  the  two  layers  must  be 
satisfied.  These  are  the  continuity  of  stress 
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where  a 


yi 


and 
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represent  stresses  acting  on  the  bottom  surface 


of  the  top  layer  of  the  beam  and  a 


and 


are  stresses  acting  on 


y2  xy2 

the  top  surface  of  the  bottom  layer  of  the  beam.  The  continuity  of 
displacements  across  the  line  of  division  consists  of 


U1  =  U2’  W1  =  w2  on  y  =  0 

Applying  boundary  conditions  (7)  and  (8)  yields 


(9) 
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From  equations  (1)  and  (3)  we  obtain 
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Integrating  the  two  equations  of  (14)  and  expressing  x  as  a  Fourier 
cosine  series 
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we  obtain 


00 

ui  =  Z  ^l(1+vPAni  +  2BnJcosh  Yy  +  (1+v  )B  Yy  sinh  yy 

1  £r,  1  Ill  Ill  1  Ill 

i  n-i 


+  I(l+v.)C  .  +  2D  ,]sinh  yy  +  (l+v.)D  .yy  cosh  yy 
x  ni  ni  l  ni 

4ta.E . 

+  - cos  Y*  +  u°  (n  =  1,3,5,...)  (16) 


-5- 


w.  =  —  2  y{(l+v.)A  .sinh  yy  +  B  .[(1  +  v.)yy  cosh  yy  -  (1-v.)  sinh  yy] 

1  E.  ,  i  ni  '  ni  l  J  i 

l  n=l 

+  (l+v^)Cn^cosh  yy  +  Dn^[(l+V^)yy  sinh  yy  “  (l~v^)cosh  yy]}sin  yx 
+  c^ty  +  w°  (n  =  1,3,5,...)  (17) 


Taking  the  point  x  =  0  and  y  =  0  as  rigidly  fixed  so  that  u^ 
#.  =  0  at  that  point,  from  equations  (16)  and  (17)  we  find 
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Applying  the  continuity  conditions  (9)  u^  =  u^  and  =  w^  along  the 
line  of  division  y  =  0  yields 
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Solving  the  eight  linear  algebraic  equations  (10)  to  (13),  (19),  and 

(20),  the  eight  unknown  coefficients  A  B  C  and  D  .  can  be 

ni  ni  ni  ni 

expressed  in  closed  forms  in  terms  of  elastic  moduli,  thermal  linear 
strain,  and  beam  dimensions.  For  numerical  solutions  these  eight 
linear  algebraic  equations  can  be  easily  solved  in  each  particular  case 
and,  once  the  eight  coefficients  are  determined,  the  following  stresses 
can  be  obtained. 
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or  .  =  -  I  y  (A  .cosh  yy  +  B  .  yy  smh  yy  +  C  .smh  yy 

yi  ,  _  *  ni  ”  ni  m 

n=l ,3, . . . 

+  D^yy  cosh  yy)sin  yx 


T  .  =  -  1  y  [A  .sinh  yy  +  B  .(sinh  yy  +  yy  cosh  yy) 

xyi  n=l ,  3 , .  . .  ni  nX 


+  C^cosh  yy  +  D^^Ccosh  yy  +  yy  sinh  yy)]cos  yx 


(22) 


00 

2 

a  .=  1  y  [A  .cosh  yy  +  B  . (2cosh  yy  +  yy  sinh  yy) 

xx  ,  „  ni  ni 

n=l ,3, . . . 

+  C^sinh  yy  +  Dni(2sinh  yy  +  yy  cosh  yy)]sin  yx 


From  equation  (22)  it  is  seen  that  one  of  the  traction-free  end  conditions 

a  .  =  0  at  x  =  0  and  x  =  £  is  already  satisfied.  The  second 
xi 

condition,  x  .  =  0  at  both  ends  of  the  beam,  cannot  be  satisfied 
xyi 

exactly  without  superimposing  additional  solutions.  However,  this 
condition  is  satisfied  in  the  sense  of  Saint  Venant's  principle.  With 
this  principle  we  may  replace  the  second  condition  by  its  statically 
equivalent  condition: 


at  x  =  0  and  x  =  £ 


(23) 


To  show  that  this  condition  is  satisfied,  we  proceed  as  follows.  Since 
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respectively.  But  the  conditions  (24)  imply  zero  resultant  shearing 
force  on  any  section,  x  =  constant,  of  the  beam.  This  follows  from 
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Thus,  the  load  on  each  end  of  the  laminated  beam  is  self-equilibrating. 
The  problem  of  exact  satisfaction  of  end  conditions  of  an  elastic  strip 
has  been  treated  by  authors  [11,12]  in  the  literature.  It  is  beyond 
the  scope  of  this  article,  however,  to  modify  the  solution  by  satisfying 
more  precisely  the  free  end  conditions  of  the  beam. 


LAMINATED  BEAMS  OF  ORTHOTROPIC  MATERIALS 


In  the  case  of  plane  stress  distribution,  the  stress-strain  and  strain- 
displacement  relations  for  orthotropic  materials  are  [10,13] 
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where 


E  ,  E  are  the  Young's  moduli  along  the  principal  directions  x 
x  y 


y;  G  =  G  ,  the  shear  modulus  which  characterizes  the  change  of  angles 
xy 


between  principal  directions  x  and  y,  and 


yx’ 


the  Poisson's  ratio. 


and 
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This  ratio  characterizes  the  decrease  in  di.cction  x  due  to  tension 

in  direction  y  with  a  similar  meaning  for  the  expression  V  ,  related 

xy 

to  V  by 
yx 
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y  xy  y  yx 


(27) 


u,  w  are  the  displacements  along  the  x  and  y  directions,  respectively, 

and  a  ,  a  represent  the  coefficients  of  thermal  expansion  in  principal 
x  y 

directions  x  and  y.  Solving  equations  (26)  for  stresses  yields 


a  =  -  [ £  +  v  e  -  (a  +  v  a  )T(x)] 

x  1-v  V  x  yx  y  x  yx  y 
xy  yx 


EE  E 

a  =  - ^ -  [ £  -t  v  £  -  (a  +  =~  v  a  )T(x)] 

y  1-v  v  1  y  E  yx  X  y  E  yx  x' 
y  xy  yx  3  y  ’  y 


The  equations  of  equilibrium  are  identically  satisfied  by  introducing 
Airy  stress  function  <j)  as 
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Txy  3x9y 


(29) 


The  substitution  of  equations  (26)  and  (29)  into  the  compatibility 
relation  yields  the  following  governing  differential  equation  for 
orthotropic  materials 
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In  the  general  case  a  complete  Fourier  series  expansion  for  T^(x)  is 
required.  However,  the  method  of  solution  remains  the  same  whether 


T\(x)  is  expanded  in  a  complete  Fourier  Series  or  a  sine  or  a  cosine 

series.  Thus,  in  the  following,  only  a  sine-series  expansion  is 
considered,  i.e., 
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Corresponding  to  the  expansion  (31),  the  particular  and  homogeneous 
solutions  for  the  ith  layer,  which  satisfy  equation  (30),  may  be 
expressed,  respectively,  as 
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in  which  A  . ,  B  . ,  C  . ,  D  .  are  arbitrary  constants  and 
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(i_  - 

'•G.  E 
1  xi 


The  general  integral  of  equation  (30)  is  the  sum  of  the  solutions  (33) 
and  (34) ,  i . e . , 


"hi 


(36) 


In  the  special  case  where  the  two  roots  \ and  of  equation  (35) 

are  real  and  equal  (i.e.,  =  1),  the  following  solution  of  equation  (30) 

should  be  used  instead  of  the  solution  (34) . 


00 

‘t’i  =  1  lAni  +  Bniy)cosh  nXty  +  (CRi  +  Dnjy)sinh  nA^Jsin  — ^  (37) 

n=l 


Having  obtained  solution  (36)  we  will  be  able  to  show  that  all  the 
boundary  conditions  (7),  (8),  and  (9)  can  be  satisfied.  Applying  the 
boundary  conditions  (7)  and  (8)  yields,  respectively. 


A  .cosh  nA.d.  +  B  .sinh  nX.d.  +  C  cosh  nu  d 
m  1  1  m  1  1  ni  i  i 


+  D  sinh  nu.d. 
ni  ri  1 


E  .a  . t  . 
yi  yi  ni 


(n  =  1,2,. ..) 


(38) 


A. (A  .sinh  nA.d.  +  B  .cosh  hX.d  ) 
1  m  11  ni  11 


+  |J.(C  sinh  np.d.  +  D  .cosh  nu.d.)  =  0, 
1  ni  1  1  ni  ri  1  ’ 


(n  =  1,2,...)  (39) 
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A  .  +  c„i  -  A  1  -  c  o  =  (£  -  E  ,)(^z)2  ,  (n  =  1 

nl  nl  n2  n2  y2  y2  n2  yl  yl  nl  nn 


A2Bn2  +  M2°n2  XlBnl  +  MlDnl’  (n  ” 


From  equations  (26)  and  (29)  we  obtain 


9ui  1 


3w.  , 

l  1 


aY 

a24>. 

1 

1 

-  V  .  — — 

ay2 

xyx  ax2 

o 

3  <J> . 
i 

a  a. 

i 

-  V  .  — — 

3x2 

yxl  ay2 

)  +  a  .T. , 


)  +  a  .T. 


Integrating  the  two  equations  of  (42)  yields 


00  2 

u.  =  —= —  1  n[(A2  +  V  .  ^r)(A  .cosh  nA.y  +  B  .sinh  nA.y) 

1  “xi  n=l  1  *»*  l2  nl  1  "> 


2  71 

+  (p.  +  v  .  — ) (C  .cosh  nu.y  +  D  .sinh  nu.y) 
i  xyi  g2  ni  ri7  ni  rr 


(E  .a  .  +  v  .E  .a  .)]cos  +  u? 
2  xi  xi  xyi  yi  yi  £  l 


oo  2 

—  1  1  2  a 

w.  =  = —  I  n{r—  (v  .A.  +  — -) (A  .sinh  nA.y  +  B  cosh  nA.y) 
l  E  .  ,  A.  yxi  i  .2  '  ni  i7  ni  i7 

yx  n=l  i  7  £ 


+  M?  (Wi  +  ^2)(Cnisinh  nMiy  +  DniCOsh  n»1iy)i,in  T 


,2,...) 

(40) 

(41) 


(42) 


(43) 


(44) 
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in  which  the  constant  u?  can  be  determined  by  taking  the  point  x  =  0 
and  y  =  0  as  rigidly  fixed  in  order  to  have  u  =  0,  =  0  at  x  =  0, 

y  =  0-  Applying  the  continuity  conditions  (9)  and 

along  the  axis  y  =  0  yields 


M  2  M  2 

iHK  +  Vxy2>Au2  +  “lb  +  Vxy2]Cn2 


Ex2  Al£  2 
xl 


*  V'*.l  *  1(^)2  *  ',xyliCnl  +  4>2  'rj  t»2<E,2»«2  *  'VW 


t  ,(E  ,a  ,  +  v  ,E  ,a  .)]}  =  0 
nl  xl  xl  xyl  yl  yl 


(45) 


X2I(A^)2  +  Vyx2lBn2  +  M2[(p2£)  +  Vyx2]Dn2  = 


+  VyxllBnl  +  M!l(ii7)2  +  VyxllDnl* 


(46) 


and  also  u°  =  u®.  (47) 

The  eight  arbitrary  coefficients  A  B  C^,  and  are  determined 

from  the  eight  linear  algebraic  equations  (38)  to  (41),  (45),  and  (46) 
and  then  the  stresses  can  be  obtained  from  the  following  expressions: 


a  . 
yi 


1  {(tt^)2  (A  .cosh  nA.y  +  B  .sinh  nA.y 
. 1  £  ni  i7  ni  x 

n=l 


+  C  .cosh  nu.y  +  D  .sinh  ny.y] 
ni  'i7  ni  ri 


+  E  .Of  .t  .  }sin 
yx  yi  ni 


nnx 

& 


(48) 


xyi 


» 

■-i  * 

n=l 


‘[A. (A  .sinh  nA.y  +  B  .cosh  nA.y) 
1  i  ni  i  ni  x 
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(49) 


+  M . (C  . 
1  ni 


.sinh  nu.y  +  D  .cosh  nu.y)]cos 
X  l'  ni  *1  £ 


2  2 

0  .  =  I  n  [A. (A  .cosh  nA.y  +  B  .sinh  nA.y) 
xi  ,  i  nx  i  ni  i 

n=l 


u2(C  .cosh  nu.y  +  D  .sinh  nu.y)]sin 
rl  ni  nx  ri  £ 


NUMERICAL  EXAMPLE  AND  CONCLUSIONS 


CASE  A:  TWO  LAYERED  BEAMS  HAVING  RIGID  BOND 


The  following  material  constants  and  beam  dimensions  are  used: 


£  =90  cm.,  hj  =  h2  =  5  cm., 

va  =  0.27  E2  =  20.69xl06  N/cm2., 

v2  =  0.33,  E2  =  6.89xl06  N/cm2., 


to  =  400°  C  =  752°  F 


a  =  6.5x10  /F° , 


a2  =  13x10  /F° , 


Figures  1  and  2  show  the  distributions  of  interlaminar  normal  and  shear 
stresses,  respectively,  along  the  longitudinal  axis  of  the  beam.  Both 
indicate  that  high  stress  intensity  occurs  in  the  end  zones  of  the  beam 
and  that  both  stresses  decay  rapidly  with  increasing  distance  from  the 
ends.  The  distribution  of  the  stresses  in  the  end  zones  within  a  distance 
equal  to  the  thickness  of  each  beam,  may  not  be  very  accurate,  especially 
for  the  shear  stress.  However,  the  results  in  figures  1  and  2  show  that 
both  stresses  do  increase  toward  the  ends  of  the  beam,  starting  from  a 
distance  greater  than  the  thickness  of  each  layer.  Thus,  it  may  be  con¬ 
cluded  that  delamination  failure,  when  it  occurs,  will  start  at  the  ends 
of  the  beam. 


Figures  3  and  4  show  the  distribution  of  axial  stress  along  the  bonding 
surface  of  the  upper  layer  and  the  lower  layer,  respectively.  The 
bonding  surface  of  the  upper  layer  is  under  tension  and  that  of  the 
lower  layer  under  compression.  Within  a  short  distance  from  the  free 
ends  of  the  beam,  the  axial  stress  of  the  upper  layer  (Fig.  3)  reaches 
its  maximum  value  and  the  axial  stress  of  the  lower  layer  (Fig.  4) 


reaches  a  value  slightly  less  than  its  maximum  value.  Both  then 
remain  constant  for  the  remainder  of  the  beam.  Other  observations 
which  could  be  made  for  Case  A  will  be  described  later  as  observa¬ 
tions  for  Case  B. 


CASE  B:  TWO  LAYERED  BEAMS  HAVING  AN  ELASTIC  BOND  BETWEEN  THE 
LAYERS  (FIG.  5) 


The  material  constants  and  beam  dimensions  for  the  two  layers  are  the 
same  as  those  used  in  the  Case  A.  The  following  constants  are  taken 
for  the  thin  bonding  material  between  the  layers: 

h3  =  .33  cm.,  v3  =  .33,  c<3  =  2.5xlO_6/F°, 

Three  cases  of  different  modulus  of  elasticity  of  the  bonding  mate¬ 
rial  are  considered,  i.e., 

(1)  E3  =  27.58xl06  N/cm?,  (2)  E3  =  4xl06  N/cm? 

(3)  E3  =  106  N/cm? 

We  note  that  the  constant  c  appearing  in  figures  1-4  and  the  con¬ 
stant  C  in  figures  6  through  13  are  related  as 


C  =  E^o^-orpc  =  44.785c. 


The  observations  and  conclusions  made  in  the  previous  Case  A  also  hold 
true  for  the  present  case.  In  addition,  other  observations  can  be  made 
from  figures  6  through  13.  Figures  6  and  8  illustrate  that  the  modulus 
of  elasticity  of  the  bonding  material  (adhesive)  affects  the  interlaminar 
normal  stress  only  in  the  narrow  end  regions,  these  regions  being  of  the 
order  of  the  thickness  of  the  beam.  A  larger  modulus  of  elasticity  in 
the  bonding  material  yields  significantly  larger  interlaminar  normal 
stress  in  the  end  regions.  However,  the  interlaminar  normal  stress  is 
virtually  unaffected  by  the  modulus  of  elasticity  of  the  adhesive  for 
the  remainder  of  the  beam.  Thus,  if  delamination  failure  should  occur, 
it  will  occur  at  the  ends  of  the  beam.  To  prevent  such  failure,  adhe¬ 
sives  having  smaller  modulus  of  elasticity  should  be  preferred. 

It  is  of  interest  to  note  that  the  modulus  of  elasticity  of  the  bonding 
material  has  only  a  slight  effect  on  the  distribution  of  the  interlaminar 
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shear  stress  as  seen  from  figures  7  and  9,  except  perhaps  in  a  very  short 
distance  from  the  ends  of  the  beam.  However,  in  this  region  the  distri¬ 
bution  of  the  shear  stress  is  not  accurate,  as  previously  pointed  out. 

An  examination  of  figures  3,  4,  10  and  11  indicates  that  the  modulus  of 
elasticity  of  the  bonding  material  has  only  slight  effect  on  the  axial 
stress  in  the  upper  and  lower  layers  of  the  beam.  On  the  other  hand, 
the  modulus  of  elasticity  of  the  bonding  material  has  significant  effect 
on  the  axial  stress  in  the  bonding  material  as  seen  from  figures  12  and 
13;  a  greater  modulus  yields  significantly  larger  axial  stress.  We 
further  note  that  the  maximum  axial  stress  in  the  upper  and  lower 
layers  of  the  beam  is  greater  than  that  of  the  interlaminar  normal 
stress.  However,  delamination  failure  of  the  bond  caused  by  the 
interlaminar  normal  stress  may  still  occur,  as  is  often  the  case  in 
practice,  if  the  bond  cannot  sustain  the  required  interlaminar  normal 
stress . 

The  computer  program  developed  for  the  numerical  calculations  is 
presented  in  the  Appendix. 


NORMAL  STRESS 


SHEAR  STRESS 


4Ett(at-a1) 


.00  .02  .04  .06  .00  .10  *w".000  .010  .020  .030  .040  -0! 

x/f  x/| 

M  146  336 

M  140  33/ 

Figure  1.  -  Interlaminar  normal  stress 
distribution. 

Figure  2.  -  Interlaminar  shear  stress 
distribution. 

AXIAL  STRESS 

AXIAL  STRESS 

4Ett(a,-o,) 

E,  *  6  89  « 10*  N/cm* 


.00  .01  .02  .03  .04  .06  .06  .07  .06 

«/f 

M  146  776 

Figure  3.  •  Axial  stress  distribution 
along  the  bonding  surface  of  the  upper 
layer  of  a  two-layered  beam. 


.00  .01  .02  .03  .04  .06  .06  .07  .06 

«/t 

M  148  776 

Figure  4.  -  Axial  stress  distribution 
along  the  bonding  surface  of  the  lower 
layer  of  a  two-layered  beam. 


Figure  S.  -  Two-layered  beam  having  an  elastic  bond  between  the  layers. 


NORMAL  STRESS 


»/t 

M  14*  7*0 

Figure  6.  -  Interlaminar  normal  stress 
distribution  along  the  bonding  surface  of 
the  lower  layer  and  the  middle  layer 
(bonding  material). 


SHEAR  STRESS 


«/« 


m  i*»  rei 

Figure  7.  •  Interlaminar  shear  stress 
distribution  along  the  bonding  surface  of 
the  lower  fayer  and  the  middle  layer 
(bonding  material). 
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NORMAL  STRESS 


SHEAR  STRESS 


»/» 


M  14®  782 

Figure  8.  -  Interlaminar  normal  stress 
distribution  along  the  bonding  surface  of 
the  upper  layer  and  the  middle  layer 
(bonding  material). 


M  148  >83 

Figure  9.  ■  Interlaminar  shear  stress 
distribution  along  the  bonding  surface  of 
the  upper  layer  and  the  middle  layer 
(bonding  material). 


AXIAL  STRESS 


M  14*7*4 

Figure  10.  •  Axial  stress  distribution 
along  the  bonding  surface  of  the  upper 
layer. 


AXIAL  STRESS 


Figure  11.  -  Axial  stress  distribution 
along  the  bonding  surface  of  the  lower 
layer. 


APPENDIX— THERMAL  STRESSES  ALONG  BEARING  SURFACES  OF  A  TWO-LAYER 
ISOTROPIC  BEAM— A  COMPUTER  PROGRAM 

#GSP,P 

PLOTTER  PEN/LIQ 
#FOR» IS  MAIN 

IMPLICIT  DOUBLE  PRECISION  (A-H,0-Z) 

C0MM0N/BLK2/  N 
COMMON/BLK3/  NUMLAY 
COMMON/BLK7/NTERMS 
COMMON/BLK12/X(301 ),Y(301) 

C0MM0N/BLK13/CC(8, 1 ) 

COKB0N/BLK20/-  A  A  (  B  ,  8  )  ,  BB  (8 , 1 ) 

CALL  MTAMDF  (AA,8,8,  '  D • ,8,8, 'GEN*) 

CALL  MTAMDF  (BB,8,1, 'D',8,1, 'GEN') 

CALL  MTADEF  (CC,8,1,  '  D '  ) 

C 

C  THIS  PROGRAM  CALCULATES  THE  STRESSES  AT  THE  INTERFACE  OF  THE  TWO 
C  ISOTROPIC  LAYERS  (YC  =  0.).  THE  NORMAL  ,  SHEAR,  AND  AXIAL  STRESSES 
C  ARE  PLOTTED  FOR  THE  UNIFORM  TEMPERATURE  C'SE. 

C 

C  READ  THE  LAST  FOURIER  TERM  TO  BE  USED.  (SHOULD  BE  ODD)  (I<l) 

READ  (5,201)  NTERMS 
CALL  INPUT 

WRITE  (6,290)  NTERMS 
I  s  1 
YC  =  .0 

DO  210  N  =  1, NTERMS, 2 
CALL  DETMAT 
CALL  SOLVE 
CALL  COEF ( I , YC ) 

210  CONTINUE 

NP01NT  =  51 
XST  ART  =  0. 

XEND  =  .1 

CALL  STRESY  ( 1  , 1 ,  NPO I  NT , YC , X S T A RT , X  END  ) 

CALL  GRAPH  (X,+2,Y,+2,NP01NT,  'NONE', 'SOLID',  '  X/L  SS','  SS* 

1, 'NORMAL', 'NORMAL  STRESSSS ', 'FULL ',' BIND  * ) 

XEND  =  .05 

CALL  STRESY  ( 2 , I , NPO I  NT , YC , XST ART , X  END ) 

CALL  GRAPH  (X,*2,Y,+2,NP01NT, 'NONE', 'SOLID', '  X/L  S$',»  $$' 

1, 'NORMAL', 'SHEAR  S TRESS $J ', 'FULL ',' BIND ' ) 

XEND  =  .08 

CALL  STRESY  ( 3 , 1 , NPOI NT , YC , XST AR T , XEND ) 

CALL  GRAPH  ( X , ♦?, Y , ♦ 2, NPOI NT ,  ' NONE  '  ,  '  SOL  I D ' , '  X/L  SS','  SS' 

1,  'NORMAL', 'AXIAL  STRFSSSS',  'FULL', 'BIND') 

CALL  STRESY  ( 3 , 2 , NPO I  NT , YC , xs T AR T , XE ND ) 

CALL  GRAPH  (X,*2,Y,+2,NP0INT, 'NONE', 'SOLID', '  X/L  S$',»  SS' 

1 , 'NORMAL  ',  'AXIAL  S T RL SS l S ' , ' F UL L ' ,  '  B I ND  '  ) 

299  STOP 

201  FORMAT  (I«) 

290  FORMAT  (  1 X,  //,  30X  ,  'THE  LAST  FOURIER  TERM  =  •  ,  I  <1 ) 

END 
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*FOR,IS  INPUT 

SUBROUTINE  INPUT 

IMPLICIT  DOUBLE  PRECISION  (A-H,0-Z) 

C  THIS  ROUTINE  READS  THE  MATERIAL  PROPERTIES  FOR  THE  2  OR  3  ISOTROPIC 

C  LAYERS  WHICH  ARE  UNDER  A  UNIFORM  TEMPERATURE  DISTRIBUTION. 

C 

COMMON/8LK1/E1 ,E2,E3,XNU1 , XNU2 , XNU3 , ALPH A  1 , ALPH A  2 , ALPH A  3 , BE T A  1 , 
18ETA2,b£TA3 
COMMON/BLK3/  NUMLAY 
C0MM0N/RLK8/XLEN 
COMMON/BLKl 0 /DELTA (3) 

C 

C  IF  ONLY  TWO  LAYERS  ARE  TO  BE  CONSIDERED,  THE  CONSTANTS  FOR  LAYER 

C  NUMBER  3  SHOULD  BE  SET  EQUAL  TO  ZERO. 

C  READ  THE  NUMBER  OF  LAYERS.  (FIRST  COLUMN  IN  INPUT  CARD) 

READ  (5,100)  NUMLAY 
WRITE  (6,180) 

WRITE  (6,190)  NUMLAY 
C 

C  READ  THE  MATERIAL  PROPERTIES.  (3F20.8) 

WRITE  (6,181) 

READ  ( 5 , ) 0  1  )  E1,E?,E3 
WRITE  (6,191)  E1,E2,E3 
READ  (5,101)  XNU1 , XNU2, XNU3 
WRITE  (6,19?)  XNU1 ,XNU2, XNU3 
READ  (5,101)  ALPHA!, ALPHA2,ALPHA3 
WRITE  (6,193)  ALPHA  1 , ALPHA?, ALPHA3 
C 

C  READ  THE  THICKNESS  OF  EACH  LAYER. 

READ  (5,101)  (DELTA(I),I=1,3) 

WRITE  ( 6 , 190)  (DELTA (I) , 1  =  1 , 3) 

C  XI.EN  s  '  THE  LENGTH  OF  THE  9EAM/PI. 

READ  (5,101)  XLEN 
WRITE  (6,195)  XLEN 
BE  T  A  1  =  DELTA (1 )/XLEN 
BETA2  =  DELT  A (2) /XLEN 
BE  T  A  3  =  DELT  A (3) /XLEN 
RETURN 

100  FORMAT  (Jl) 

101  FORMAT  (3E20.8) 

180  FORMAT  (  •  1  • , 37X THERMAL  STRESSES  IN  LAMINATED  BEAMS  OF  ISOTROPIC 
1MATERIALS* ,//) 

181  FORMAT  (IX,//, 05X,  'LAYER  1  * , 1 3X , »L AYER  2  1 ,  1  3X ,  '  L  A  YER  3') 

190  FORMAT  (JX,//,«5X,  'NUMBfcR  OF  ISOTROPIC  LAYERS  =  ',11,//) 

191  FORMAT  (IX,//, 10X,  'YOUNG  S  MODULOUS • , 1 2 X , E 1 5 . 5 , 2 ( 5X , E 1 5 . 5 ) ) 

192  FORMAT  (IX,//, 10X, 'POISSON  S  R  A  T  I  0  • ,  1  3  X  ,  E  1  5 . 5 , 2  ( 5  X  ,  F.  1  5 . 5  )  ) 

193  FORMAT  (IX,//, 10X,  'THERMAL  EXPANSION  COEF  '  ,6X,E 15.5*2(5X,E 15.5) ) 

190  FORMAT  (IX,//, 10X, 'THICKNESS  8X , E 1 5 . 5 , 2 (5X , E 1 5 . 5 ) ) 

195  FORMAT  (1X,////,30X, ’ (BE  AM  LENGTH)  /  PI  =  ',E15.5) 

END 
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*FOR,IS  DETMAT 

SUBROUTINE  DETMAT 

IMPLICIT  DOUBLE  PRECISION  (A-H,0-Z) 

C  THIS  SUBROUTINE  DETERMINES  THE  MATRICIES  FOR  THE  TWO  LAYER  ISOTROPIC 
C  AND  UNIFORM  TEMPERATURE  CASE 
C 

COMMON/BLK 1/E  1 , E2,E3, XNUl , XNU2, XNU3, ALPHA1 , ALPHA2, ALPHA3,BETA1 , 

1  BET  A2, BET  A3 
COMMON/BLK2/  N 
C0MK0N/BLM/A(8,8),B(8) 

C 

IF  (N.GT.l)  GO  TO  10 
D 1  =  E2/E1 
D2  =  1.  +  XNUl 
D3  =1.  +  XNU2 
DO  1  I  =  1,8 
H ( I )  =  0. 

DO  2  J  =  1,8 
A ( I , J )  =  0. 

2  CONTINUE 
1  CONTINUE 
A  (  1  , 1 )  =  1. 

A ( 2 , 5 )  =  1. 

A  (  3 , 3  )  =  1  . 

A («,7)  =  -1  . 

A(5,l)  =  1  . 

A ( 5  »  5)  =  -1. 

A(fe,3)  =  1  . 

A(b,<J)  =  1  . 

A ( 6 , 7  )  =  -1  . 

A  (  6 , 8 )  =  -1 . 

A (8, 3)  =  01 *D2 
A(7,l)  =  A ( 8 , 3  ) 

A  (  7 , 2  )  =  2  .  *  D  1 
A (7,5)  =  -1 .*03 
A  (  7  ,  b )  =  -2. 

A(8,<l)  =  -1  ,*D1*(1  .  -  XNUl  ) 

A ( 8 , 7 )  =  A (7,5) 

A ( 8 , 8 )  si.-  XNU2 
10  XN  =  N 

BIN  =  XN*BETA1 
B2N  =  XN*BFTA2 
TH1  =0TANH(B1N) 

TH2  =DTANH(B2n) 

A  ( 1 , 2  )  =  B  1  N  *  T  H  1 
A  ( 1 , 3 )  =  TH1 
A  (  1  ,  <1 )  =  BIN 
A  (  3 , J )  =  TH1 
A ( 3 , 2 )  =  BIN  ♦  TH1 
A  (  3,  <1 )  -  1.  ♦  A  (1 ,2) 

A ( 2, 6 )  =  B2N*TH2 
A  f  2, 7 )  =  -1 . *  T  H2 
A (2,8)  =  - 1 . *  0  2  N 
A  (  u , 5 )  =  TH2 
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A  ( 4 , 6  )  r  B2N  t  TH2 

A(4,8)  =  -1.  -  (B2N*TH2) 

B ( 7 )  =  l./XN 

RETURN 

END 


•FOR, IS  SOLVE 

SUBROUTINE  SOLVE 

IMPLICIT  DOUBLE  PRECISION  (A-H,0-Z) 

C  THIS  SUBROUTINE  SOLVES  THE  MATRIX  EQUATIONS  FOR  A  GIVEN  N  FOR  THE 
C  TWO  OR  THREE  LAYER  CASE.  THE  REQUIRED  INPUT  IS  THE  FOURIER  NUMBER  N » 

C  THE  NUMBER  OF  LAYERS  NUMLAY,  AND  THE  MATRIX  A  AND  VECTOR  B  FROM 

C  SUBROUTINE  DETMAT.  THE  MACC  SUBROUTINE  MTSOLV  IS  UTILIZED. 

C  THE  SOLUTION  VECTOR  IS  CC(K,1). 

COMMQN/BLK2/  N 
COMMON/BLK3/  NUMLAY 
C0MM0N/BLK4/AMAT (8,8) , 8VEC (8) 

CQMM0N/RLK!3/CC(8, 1  ) 

COMMON/8LK20/  AA(8,B),BB(8, 1) 

MATSIZ  =  4*NUMLAY 
DO  410  I  =  1, MATSIZ 
8B(I,1)  =  BVEC(I) 

DO  420  J  =  1, MATSIZ 
AA(I,J)  =  AMAT(IrJ) 

420  CONTINUE 
410  CONTINUE 

CALL  MTSLVD  (AA,BB,CC, IRET) 

299  RETURN 
END 


“FOR, IS  COEF 

SUBROUTINE  COEF(J,YC) 

IMPLICIT  DOUBLE  PRECISION  (A-H,0-Z) 

C  THIS  SUBROUTINE  DETERMINES  THE  FOURIER  COEFFICIENTS  FOR  THE  STRESSES 
C  IN  THE  I  TH  LAYER  ON  Y  =  0 . 

C0MM0N/8LK2/  N 
COMMON/8LK8/XLEN 
COMMON /BLK 1 3/BB  (8, 1  ) 

COMMON/BLK 1 4/COEFN(  1800) ,COEFS( 1800) , COEF AC  1800,2) 

ZN  =  N 

M  =  (N  ♦  1 )/2 
COEFN(M)  =  BB  ( 1  ,  1  ) 

COEF  S (M)  =  B  B ( 3  ,  1  )  +  eB(4, 1  ) 

COEFA  ( M , 1 )  =  BB ( 1 , 1 )  ♦  2.*BB(2,1) 

COEFA  ( M , 2 )  =  RB(5,J)  +  2.*BB(6,1) 

RETURN 

END 


FFOR/IS  STRESY 

SUBROUTINE  STRESY  ( I T YPF , I , NPO I  NT , YC , XS T AR T , XE NO ) 
IMPLICIT  DOUBLE  PRECISION  (A-H,0-Z) 

-24- 


o  n 


1 


C  THIS  SUBROUTINE  COMPUTES  THE  NORMAL,  SHEAR,  OR  AXIAL  STRESSES  ON 
C  THE  LINE  Y  =  0,  AT  NPOINT  NUMBER  OF  POINTS  FROM  X  =  XST  ART  TO  X  = 

C  XEND.  (X  :  0  TO  1  IS  THE  ENTIRE  LENGTH  OF  THE  BEAM)  IF  I  TYPE  =  1, 

C  THE  NORMAL  STRESS  IS  COMPUTED,  IF  ITYPE  =  2,  THE  SHEAR  STRESS  IS 

C  COMPUTED,  OR  IF  ITYPE  =  3,  THE  AXIAL  STRESS  IS  COMPUTED.  THE  (X,Y) 

C  COORDINATES  ARE  RETURNED  IN  THE  VECTORS  X  AND  Y.  THESE  STRESSES 
C  APE  CALCULATED  FOR  THE  I  TH  LAYER. 

COMMON/BLK7/NTERMS 
COMMON/BLK6/XLEN 
COMMON/8LK12/X(301 ) , Y(301 ) 

C0MM0N/BLK14/C0EFN( 1800) ,COEFS ( 1800) ,COEFA (1800,2) 

IF  (ITYPE. EQ. 2)  GO  TO  600 
IF  (ITYPE. EO. 3)  GO  TO  611 
WRITE  (6,691)  YC 


GO  TO  601 

600  WRITE  (6,692) 

GO  TO  601 

611  WRITE  (6,671) 

601  WRITE  (6,693) 

WRITE  (6,694) 

THIS  IS  DONE  SO 
=  1,  THIS  ROUTINE  WILL 
XPOINT  =  NPOINT  -  1 


YC 

YC 

I , NPOINT 
XSTART,XEND 
BOTH  END  POINTS 

fail. 


WILL  BE  PLOTTED  .  HOWEVER,  IF  NPOINT 


720 


700 


756 


PI  =  3.1415926536 
IF  (ITYPE. EQ. 2)  GO  TO  755 
IF  (ITYPE. EO. 3)  GO  TO  756 
DO  700  J  =  1, NPOINT 
ZJ  =  J  -  1 

X(J)  =  ZJ*(XEND  -  XSTART)/XPOINT  +  XST  ART 
ARG  =  X(J)*PI 
XSUM  =  0. 

DO  720  N  =  1 ,NTERMS,2 
M  =  (N  +  l)/2 
XN  =  N 

XSUM  =  XSUM  +  COEFN(M) * DS I N ( ARG * XN ) 
CONTINUE 
Y(J)  =  -l.*XSUM 
CONTINUE 
GO  TO  770 

DO  701  J  =  1, NPOINT 
ZJ  =  J  -  1 


ZJ  3  J  -  1 

X(J)  3  ZJ*(X£NO  -  XSTART  J/XPOINT  t  XSTART 
ARC  =  X(J)*PI 
XSUM  s  0. 

00  730  N  =  1,NTERMS,2 
M  s  (N  ♦  l)/2 
XN  =  N 

XSUM  s  XSUM  +  C0EFS(M)*DC0S(ARG*XN) 

730  CONTINUE 

Y ( J )  s  -l.*XSUM 
710  CONTINUE 
770  RETURN 

69<l  FORMAT  (1 X, ///, 20X, > THIS  STRESS  IS  TO  BE  PLOTTED  FROM  X  s  *,F6.3,' 
1  TO  SF6.3) 

691  FORMAT  (IX, /////, 20X,  '  THE  NORMAL  STRESS  ON  THE  LINE  Y  s  ',F6.3, 

1  •  IS  PLOTTED  • ) 

692  FORMAT  (IX, /////, 20X, 'THE  SHEAR  STRESS  ON  THE  LINE  Y  3  ',F6.3, 

1  •  IS  PLOTTED  ' ) 

671  FORMAT  ( 1 X ,/////, BOX T HE  AXIAL  STRESS  ON  THE  LINE  Y  3  «,F6.3, 

1  '  IS  PLOTTED  ' ) 

693  FORMAT  (20X,'FOR  LAYER  NUMBER  ',Ilf'  .  THE  NUMBER  OF  POINTS  PLOT 


1  TED  IS  '14,'  . ' ) 

780  FORMAT  (13) 

781  FORMAT  (F10. 5) 

790  FORMAT  ( 1 X , 0 OX ,  'NUMBER  OF  POINTS  TO  BE  PLOTTED  s  »,I3) 

791  FORMAT  (1X,*10X,'  LENGTH  OF  BEAM  TO  BE  PLOTTED  s  *,F6.3) 
END 
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